Non-Kondo mechanism for resistivity minimum in spin ice conduction systems 
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We present a mechanism of resistivity minimum in conduction electron systems coupled with 
localized moments, which is distinguished from the Kondo effect. Instead of the spin-flip process 
in the Kondo effect, electrons are elastically scattered by local spin correlations which evolve in a 
particular way under geometrical frustration as decreasing temperature. This is demonstrated by the 
cellular dynamical mean-field theory for a spin-ice type Kondo lattice model on a pyrochlore lattice. 
Peculiar temperature dependences of the resistivity, specific heat, and magnetic susceptibility in the 
non-Kondo mechanism are compared with the experimental data in metallic Ir pyrochlore oxides. 

PACS numbers: 71.10.Fd, 71.20.Be, 71.20.Eh, 72.15.-v 



The interplay between conduction electrons and local- 
ized moments is one of the most fundamental subjects in 
strongly correlated electron systems. It leads to a vari- 
ety of fascinating phenomena in electronic and magnetic 
properties. Recent intensive studies have revealed that 
anomalous magneto-transport properties emerge from 
the coupling to localized magnetic moments 0, Q . 

Kondo effect is a basic concept to understand the prop- 
erties of electron-spin coupled systems. As its simplest 
manifestation, the electrical resistivity shows a minimum 
as decreasing temperature (T) in a metal with magnetic 
impurities. Kondo successfully explained the origin of 
resistivity minimum in terms of enhanced spin-flip scat- 
tering by magnetic impurities at low T Q . After his pio- 
neering work, the theory was elaborated to give a unified 
description of the system through a single energy scale, 
the Kondo temperature Tk 

Recently, resistivity minimum has drawn a renewed 
interest since it was observed in several metallic Ir py- 
rochlore oxides, such as Pr 2 Ir 2 07 @ and Nd2lr2C>7 un- 
der pressure Q. These compounds include conduction 
electrons in Ir 5d orbitals interacting with magnetic mo- 
ments from rare-earth 4/ localized electrons. The situa- 
tion can be well characterized as a Kondo lattice system, 
and importance of quantum fluctuations in 4/ moments 
was pointed out [!, 0-QiJ- Nevertheless, it is still un- 
clear whether the Kondo scenario is straightforwardly 
applied to the observed resistivity minimum. Firstly, 
the local ground state of 4/ multiplct is a non-Kramers 
doublet separated from excited states by a fairly large 
gap ~ 160 K jn|. This leads to a strong easy-axis 
anisotropy of 4/ moments along the (111) directions in 
the pyrochlore lattice (see Fig. [1} . The strong anisotropy 
considerably reduces Tk, and hence, it is unclear whether 
the resistivity minimum at high T <~ 40 K is explained 
by the Kondo effect. Secondly, the magnetic suscepti- 
bility shows a diverging behavior below 40 K, in sharp 
contrast to a saturation in the Kondo effect @. Thirdly, 
the specific heat shows a peak at ~ 2 K much lower than 
40 K suggesting the absence of entropy release due 
to the Kondo effect around 40 K. 

In this Letter, we propose an alternative mechanism for 



the unconventional resistivity minimum. Instead of the 
Kondo effect, we focus on the role of anomalous spatial 
magnetic correlation imposed by geometrical frustration. 
In these systems, it is inferred that localized moments 
develop spin-ice type local correlation [f| Q; that is, at 
four vertices in each tetrahedron in the pyrochlore lattice, 
two spins point inward and the other two point outward, 
without showing long-range ordering [12, [H|. In the 
present study, we show that conduction electrons inter- 
acting with the spin-ice type moments exhibit a resistiv- 
ity minimum at a characteristic temperature T m ; n , below 
which the spin-ice correlation develops. The magnetic 
susceptibility and the specific heat also show character- 
istic T dependences qualitatively different from those in 
Kondo impurity systems. The results are discussed in 
comparison with the experimental data for Pr2lr207 and 
Nd 2 Ir 2 7 . 

We start with a Kondo lattice model on a pyrochlore 
lattice (Fig. [1} , whose Hamiltonian is given by 



* = -* £ (4 



H.c. 



s., (1) 



Here, the first term describes the kinetic energy of elec- 
trons. The sum is taken over the nearest-neighbor 
(n.n.) sites. The second term represents the coupling be- 
tween conduction electron spins and localized moments 
Si (<x is the Pauli matrix). We consider the limit of 
strong easy-axis anisotropy along the (111) directions for 
the localized moments, namely, Si is a classical Ising spin 
(| Si | = 1) parallel to the local easy axis at each four- 
sublattice site, as shown in Fig.[IJa). In this limit, quan- 
tum fluctuation of the charge-spin coupling is completely 
suppressed: the sign of J does not matter, and the pos- 
sibility of Kondo effect is excluded. Hereafter, we take 
t = 1 as an energy unit, and the lattice constant of the 
16-site cubic unit cell as a length unit . 

To study this model, we adopt the cellular dynamical 
mean-field theory (CDMFT) [lj]. Through the previous 
studies of a double-exchange model, DMFT was shown 
to be a powerful tool for describing the physics of con- 
duction electrons coupled to thermally fluctuating clas- 
sical spins [151 ] . To take into account the peculiar spatial 
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FIG. 1. (color online), (a) A cubic unit cell of a pyrochlore 
lattice. The left panel shows the four sublattices with the 
(111) easy axes. The right panels show 4- and 16-site clusters 
used in the CDMFT calculations, (b) and (c) illustrate the 
magnetic ordering patterns for two-in two-out and all-in/all- 
out ordered states, respectively. 



correlations under the geometrical frustration in the py- 
rochlore structure, we analyze the model ([T]) by using the 
CDMFT with clusters composed of N c = 4 and 16 sites 
as shown in Fig. (Ua). We start with an initial Weiss 
field (icp), where e p — (2p + l)irT is the Matsubara 
frequency, and r\ = (a, s) takes 2N C values corresponding 
to the site index in each cluster a and spin s. The local 
Green's function is obtained by 

G(ie p ) = (G^(ie p )) = Tr {Sa} G^ (te p )P({S a })(2) 

where P({S a }) = e -Se«({S«}) /Tr {Sa}e - s * ff « s <"» with 
e -s off ({s Q }) = l\ p dct[~G^(ie p )}, and (3 = 1/T. Here, 

G { *; ] (ie P ) is given by G^(ie p ) = \g~\iep) + J(<r • 
S a )(5 aQ '] _1 for each spin configuration {S Q }. (The 
hat indicates a 2N C x 2N C matrix). We calculate 
the trace in Eq. @ exactly by enumerating all 2 N ° 
spin configurations with the help of fast update algo- 
rithm (l6| . Hence, we arc free from any approxima- 
tions or statistical errors in solving the impurity clus- 
ter problem. A closed set of the self-consistent equa- 
tions of CDMFT is given by Eq. the Dyson equa- 
tion, Tj(ie p ) = g (ie p ) — G _1 (ie p ), and the relation 
between the lattice and the impurity Green's function, 
G{ie p ) = £ Ek[* e P - *( k ) + M - S(iep)]" 1 . Here t(k) is 
the Fourier transform of the hopping matrix in Eq. ([1} , 
u is the chemical potential, and N is the total number 
of clusters in the whole lattice. Dynamical quantities 
arc calculated via the retarded Green's function obtained 
simply by replacing ie p by e + iS in Eq. The electri- 
cal conductivity is calculated by the Kubo formula with 
evaluating current-current correlation functions on the 
real-frequency axis. 

Figure HJa) shows the phase diagram obtained by 
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FIG. 2. (color online), (a) Particle density (n) and charge- 
spin coupling (J) dependences of the transition temperature 
(T c ). The thin dashed lines indicate the phase boundaries 
at T = 0. The gray area shows the phase separated region. 

(b) n dependences of n.n. component of RKKY interaction 
(Jrkky), Curie-Weiss temperature (Tew), and T c at J = 2. 

(c) and (d) show T dependences of the probabilities for each 
spin configuration in a tetrahedron at n = 0.05 and 0.40, 
respectively: P22, P31, and P40 correspond to the two-in two- 
out, three-in one-out or one-in three-out, and all-in or all-out 
configuration, respectively. The results of N c = 4 and 16 are 
plotted at J = 1. The dashed curve in (d) shows the order 
parameter for the all- in/ all- out state. 



CDMFT. We focus on the density region less than half 
filling, n = J2icr( c ia c i<?) /NN C < 0.5. Hereafter, we as- 
sume a single tetrahedron as a magnetic unit cell [l7| . 
We obtain two dominant ground states in the present 
paramter range: One is a two-in two-out ordered state in 
the low-n region [Fig.QJb)], and the other is an all-in/all- 
out ordered state near half filling [Fig.QJc)]. Comparison 
of the ground-state energies between two phases gives the 
phase boundary at n ~ 0.2 with a phase-separated re- 
gion, as shown in Fig. Ufa). The gradation and contour 
lines represent the transition temperature T c for each or- 
dering calculated by CDMFT with N c = 4 JS]. The 
results are almost independent of the cluster size, as wc 
will see below. 

In order to understand the phase diagram, we compare 
the n.n. component of the Ruderman-Kittcl-Kasuya- 
Yosida (RKKY) interaction, Jrkky [III, the Curie- Weiss 
temperature Tew, and T c in Fig. &b) for J = 2. Here, 
Jrkky is calculated from the eigenstatcs of Hamiltonian 
([T]) at J = 0, by the second-order perturbation in terms 
of J/t. Tqw is obtained by the Curie- Weiss fit for the 
magnetic susceptibility as x ^ (T — Tqw ) i n the range 
of 0.5 < T < 1.0 [see also Fig. G|d)] 0. The ground 
state is well understood from Jrkky- Jrkky is ferro- 
magnetic in the low-n region, reflecting a small Fermi sur- 
face. The ferromagnetic interaction between neighboring 
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localized moments favors the two-in two-out configura- 
tion. As n is increased, however, Jrkky changes its sign 
at n ~ 0.2 and becomes antifcrromagnctic, which favors 
the all- in/all-out configuration appearing for n > 0.2. On 
the other hand, Tew has less correlation with the order- 
ing pattern. This may be attributed to the long-range 
and oscillating nature of RKK Y interaction [2l[ . 

Although the sign of Jrkky well accounts for the dom- 
inant phases at T = 0, the magnitude of Jrkky does not 
explain large difference of T c between two phases. T c for 
the two-in two-out ordering is largely suppressed com- 
pared to that for the all- in/all-out ordering, despite the 
magnitude of Jrkky is similar in the two regions, as 
shown in Fig. [2fb) . The contrasting T c may be a mani- 
festation of the frustration effect. Namely, ferromagnetic 
Jrkky is underconstraint and results in macroscopic de- 
generacy among different configurations of two-in two- 
out tetrahedra, as is well known in spin ice [HI, [l3| : A 
true long-range ordering is driven by further-neighbor 
RKKY interactions, which are smaller than the n.n. 
Jrkky- On the other hand, the all-in/all-out configura- 
tion is free from the frustration and a unique long-range 
order is selected by the n.n. antiferromagnetic Jrkky 
alone. Thus, the geometrical frustration plays an impor- 
tant role on the suppression of two-in two-out ordering 
in the low-n region. 

The suppression of T c leads to emergence of peculiar 
liquid-like behavior above T c . In Figs. [2jc) and HJd), 
we plot the weights of local spin configurations within a 
tetrahedron for n = 0.05 and 0.40 at J = 1, respectively. 
P22, P31, and P40 are the probabilities that we find two-in 
two-out, three- in one-out or one-in three-out, and all-in 
or all-out configurations in a tctrahedral unit, respec- 
tively. The results of N c = 4 and 16 coincide with each 
other to a remarkable extent. For n = 0.40, a transition 
with all-in/all-out ordering occurs at T c ~ 0.05, and ac- 
cordingly, P40 shows a kink-like singularity and increases 
rapidly below T c . On the other hand, for n = 0.05, 
T c is suppressed below 0.001, despite that P22 starts to 
increase considerably at T ~ 0.1 and exceeds 0.999 at 
T = 0.001. Thus, the two-in two-out spin-ice type local 
correlation is developed from much higher T than T c , in 
sharp contrast to the all-in/all-out case. The growth of 
P22 without long-range ordering is a clear signature of 
emergence of spin liquid state. 

We find that the peculiar liquid state exhibits an in- 
teresting transport property. Figure 03a) shows T de- 
pendence of the electrical resistivity p at n = 0.05 for 
three different values of J, obtained with N c = 4. The 
resistivity exhibits a clear minimum at T = T m - m and 
subsequent upturn below T m ; n . The resistivity upturn 
is somewhat counterintuitive. Usually, in electron-spin 
coupled systems, magnetic correlations develop at low T 
so as to lower the kinetic energy and restore the coherent 
motion of electrons, which leads to the reduction of resis- 
tivity. The mechanism is widely seen in a ferromagnetic 
ordering of the double-exchange systems [l5[ . Then, the 
question is what brings about the resistivity upturn be- 




FIG. 3. (color online). T dependences of (a) resistivity p, (b) 
P22, (c) specific heat C v , and (d) magnetic susceptibility \- 
The data are obtained at n = 0.05 by CDMFT with N c = 4 
cluster, p is measured in unit of h/8ne 2 , where h is Planck 
constant and e is elementary electronic charge. To facilitate 
the comparison, the data for J = 0.7 and 0.4 are multiplied 
by 1.8 and 5.0, respectively. The inset of (a) shows the com- 
parison of five different levels of approximation. In (d), the 
Curie- Weiss fitting xcw and "Kondo" susceptibility xk are 
shown for comparison. The inset shows the plot of x an d X~ 
in the linear T scale. 



low T m j n in the present pyrochlore system. 

In order to gain an insight into the mechanism, we 
compare p calculated at several different levels of approx- 
imation. The results are shown in the inset of Fig. . 
In the first place, p shows only a small change by ig- 
noring the vertex correction T, which implies that the 
self-energy effect or enhanced scattering rate plays a cen- 
tral role in the resistivity minimum. To clarify the ef- 
fects of spin correlations via the self-energy, we use the 
perturbation up to the second order in terms of J/t as 
= J 2 g n ri'(Sa • S a '). When neglecting spatial spin 
correlation in £ by taking {S a • S Q ') = 5 aa ', p decreases 
monotonically as lowering T. (This corresponds to the 
second-order Born approximation.) Alternatively, when 



4 



we use the CDMFT result of (S a ■ S a >) in £ to incorpo- 
rate the effect of growing spin ice correlation (we call this 
scheme the "correlated Born approximation"), p shows 
clear upturn, qualitatively consistent with the nonpertur- 
bative CDMFT result. This striking difference by consid- 
ering short-range spin correlations corroborates the cru- 
cial importance of spin ice correlation in the resistivity 
upturn. Indeed, T m - m roughly corresponds to the onset 
of two-in two-out correlation J22, as shown in Figs. [3ta) 
and[3{b). We also note that the result for N c = 16 gives 
a steeper upturn than that for N c = 4, suggesting that 
spatial correlations beyond the n.n. sites enhance the 
scattering of electrons. 

The above considerations clearly show that the resis- 
tivity minimum in the present system is caused by the 
peculiar local spin correlation growing under the strong 
geometrical frustration. Here, the lowering of conduc- 
tion electron energy gives rise to the enhancement of 
two-in two-out local spin correlations, leading to a for- 
mation of spatially disordered configurations of two-in 
two-out tetrahedra. Our results indicate that the "spin- 
ice manifold" acts as a rather stronger scatterer for elec- 
trons, than simple paramagnetic moments. Below T c , 
two-in two-out tetrahedra are aligned and the degener- 
acy of the manifold is lifted, leading to a sharp drop of p. 
Similarly, an application of magnetic field above T c also 
causes a reduction of p, i.e., negative magnetoresistance 
(not shown) . The coexistence of seemingly contradicting 
features — reduction of conduction electron energy and 
enhanced scattering rate — is indeed a manifestation of 
unusual physics in the itinerant frustrated system. 

Now we discuss T dependence of thermodynamic 
quantities. Figure (3Jc) shows the specific heat C v = 
TdS/dT, calculated from the spin entropy S = 
-Xi { s Q } P ({S Q })log.P({S Q }) [H. The results are ob- 
tained with N c = 16. C v does not have any significant 
feature at T ~ T m ; n , but shows a broad peak at a much 
lower T. The peak temperature roughly corresponds to 
the inflection point of P22 in Fig. [S^b), suggesting that 
the entropy is released at the saturation of spin ice corre- 
lation, not at the onset. On the other hand, the magnetic 
susceptibility x continuously increases as T decreases, as 
shown in Fig. Eld). The high-T behavior is well fitted 
by the Curie- Weiss law, xcw = C/(T — Tew) [see in- 
set of Fig. Efd)]. Tew depends on n as well as J, as 
shown in Fig. Et^b). As lowering T, x deviates from xcw 
below T m i n , while it continues to diverge weaker than 
Xcw, as shown in Fig. EJd). For comparison, we plot a 
mimic of the high-T susceptibility of the Kondo system, 
Xk = C/(T+T K ), with taking a small enough T K = 0.01. 
The diverging x as well as the broad peak of C v at low T 
is in marked contrast with the canonical Kondo behav- 



iors. 

Finally, let us discuss our results in comparison with 
the experimental data of Pr2lr207 [f| and Nd2lr207 un- 
der pressure @ . Our scenario of resistivity minimum due 
to the growing spin-ice local correlation explains many 
experimental features, at least qualitatively, such as the 
large separation of T scales between the resistivity min- 
imum and specific heat peak, and the diverging behav- 
ior of the magnetic susceptibility Q. These features are 
expected to remain robust down to a smaller realistic 
value of J/t. They will also be unaltered even when the 
exchange interaction between localized moments is in- 
cluded in the model, which is small but a finite in real 
compounds. Compared to the canonical Kondo effect, 
our scenario is free from the difficulties mentioned in the 
introduction, and therefore, potentially gives a unified 
explanation for the resistivity upturn in experiments. It 
also naturally explains the resistivity drop by long-range 
ordering at a lower T than T m j n observed in Nd2lr20y 
under pressure @. A smoking-gun experiment is a pho- 
tocmission spectroscopy; a Kondo resonance peak, which 
is naturally expected to develop below T m ; n in the canon- 
ical Kondo systems, is absent in our scenario. 

In summary, we have studied the transport properties 
of conduction electrons interacting with spin-ice type mo- 
ments. We found that the resistivity takes a minimum 
at a characteristic temperature, below which local spin- 
ice correlations begin to grow under geometrical frustra- 
tion. Our results agree with the experimental results 
of Pr2lr207 and Nd2h"207 in many respects, support- 
ing our scenario as a promising mechanism of resistiv- 
ity upturn in these compounds. The family of Ir py- 
rochlorcs exhibits metal-insulator transitions by chemi- 
cal substitutions [23|, [H| or external pressure (g, [2f| , in 
which Coulomb repulsion and spin-orbit interaction may 
play an important role. Furthermore, in P^b^Oj, an 
anomalous Hall effect without magnetic ordering was re- 
ported at very low T, and the importance of quantum 
fluctuations of Pr moments was discussed It is in- 

teresting to extend the present study for comprehensive 
understanding of these phenomena by including these el- 
ements neglected in our model. It is also intriguing to ask 
how our scenario is universal. Other types of spatial fluc- 
tuations such as orbital and lattice degrees of freedom, 
may lead to resistivity upturn as well as other fascinat- 
ing transport phenomena under severe frustration. The 
present results give a firm ground for such extensions. 
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